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INTRODUCTION 


A number of upcoming space projects - for example, Langley Research Center’s 
SAGE II and HALOE projects - make use of instruments carried aboard Earth orbital 
satellites which view the Sun and measure the intensity of solar radiation as the Sun 
"rises" or "sets" with respect to the spacecraft* The atmospheric optical path 
through which the spacecraft views the Sun varies in length during such an event and 
attenuates the measured solar radiation to varying degrees. Consequently, by the use 
of established inversion techniques, the vertical distribution of optically active 
constituents along the ray path can be mapped* 

A number of geometric parameters defined by the relative positions of the space- 
craft, the Earth, and the Sun are used both in mission planning and during mission 
operations and data reduction. Many of these parameters are highly dependent on the 
orientation of the orbital plane with the Earth-Sun direction. 

In the mission planning stages, for example, a number of these parameters are 
used to determine the maximum required excursion of instrument pointing angles - 
elevation and azimuth angles - and their rates of change for proper sizing and place- 
ment of structures such as gimbals and drive mechanisms and for imposing accuracy 
requirements on the design performance of these structures. The tangent height, 
defined as the minimum altitude above the Earth's surface of the line of sight from 
the spacecraft to some point of the solar disk, is also an important parameter. In 
addition to determining the geographic "location" of the measurement, it is used to 
define the mission event time, which is the time required to pass between two given 
tangent heights. This time determines the data collection period, which in turn can 
set requirements on data collection rates to insure that sufficient data are taken 
during the event. This time may also affect the thermal design of the instrument or 
its peripheral support components, as it determines the total time the instrument 
hardware might be exposed to the direct rays of the Sun. 

Traditionally, for elliptic orbits over an oblate Earth, these parameters are 
computed by dividing the spacecraft orbit into several hundred steps, computing these 
parameters at each step, and then interpolating in the resultant tables for the 
required values. This procedure can be rather expensive in terms of computer time, 
as over a typical 3- or 4-year mission, it may be necessary to compute thousands of 
orbits to determine the maximum design ranges of these instrument parameters and to 
predict the time and location of future events. For mission design purposes, how- 
ever, the analysis can be greatly simplified if one is willing to assume that only 
the secular perturbations produced by the oblate gravity field are important and that 
the short-period and long-period perturbations, which are generally an order of mag- 
nitude or more smaller than the secular perturbations, can be neglected. Using these 
assumptions. Brooks (1977) presents a concise but thorough review of orbit mechanics 
from a mission-planning point of view. Brooks (1980) applies these results to the 
design of solar occultation missions. This paper concentrates on a detailed examina- 
tion of the rise or set problem, which is the operational mode of many Earth satel- 
lite experiments and hence complements the material presented in the two papers by 
Brooks. The inclusion of the appendix in the present paper eliminates the need for 
reliance on external data such as almanacs or ephemerides. 
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SYMBOLS 


A 2 azimuth angle, deg 

a semimajor axis of orbit, km 

E eccentric anomaly, deg 

E^ elevation angle above or below local horizontal, deg 

e eccentricity 

e r unit vector along radius of spacecraft 

e g unit vector to some part of solar disk 

e 1 unit vector to the Sun 

e 2 unit vector in orbital plane 

e^ defined by equation (13) 

GMT Greenwich mean time 

i inclination of orbital plane to Earth equator, deg 

M mean anomaly, deg 

n mean angular rate of spacecraft in its orbit, rad/sec 

unit vector in orbital plane pointing toward periapsis 
,q,w direction cosines of e g vector along p, Q, and w, respectively 

unit vector in orbital plane pointing toward a true anomaly of 90° 
magnitude of spacecraft radius vector, km 
radius of the Earth, km 

. vector which forms a triangle with R. and p , km 

c t t 

STG sidereal time at Greenwich 

STGO Greenwich sidereal time at 0 hours GMT 
T last time of periapsis passage 

t time 

v true anomaly, deg 

w unit vector normal to plane of orbit along positive direction of angular 

momentum vector 
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A 

Xq component of radius vector along p, km 

Y Q component of radius vector along Q, km 

oi g right ascension of the Sun, deg 

a t right ascension of p t , deg 

6 angle between Sun vector and orbital plane, deg 

Y angular arc traced by subtangent point during rise or set event, deg 

S g declination of the Sun, deg 

6^ declination of p fc , deg 

e obliquity of the ecliptic, deg 

0 angle between radius vector and Sun vector, deg 

local longitude (positive east), deg 
p magnitude of projection of £ onto e 2 -e 3 plane, km 

a . A 

a x component of e g m -w direction 

A A A 

a_, component of e o in W x e direction 

Y s r 

ft longitude of ascending node, deg 

a) argument of periapsis, deg 

co g angular rate of the Earth about the Sun, rad/sec 

Subscripts : 

max maximum 

min minimum 

t sunrise or sunset event 

Superscripts : 

h hours 

m minutes 

s seconds 

A dot (•) over a symbol denotes a derivative with respect to time? a caret ( A ) 
over a symbol denotes a unit vector; and an arrow (-►) over a symbol denotes a 
dimensioned vector* 
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This paper presents an analytical method for computing the parameters pertinent 
to sunrise or sunset events (hereafter called events). No interpolation is involved, 
and only the parameters of interest are computed at the time and altitude conditions 
imposed by mission requirements. Two assumptions are made (1) the Earth is spher- 
ical, and (2) refraction is neglected. Both of these assumptions, however, can be 
removed by iterating on the solutions presented herein. A general elliptical orbit 
is assumed initially. However, once the basic equations are derived, a restriction 
to circular orbits is introduced, which greatly simplifies the resulting equations 
and permits a ready interpretation for extension to the more complex elliptic orbit 
geometry. 


ANALYSIS 

One of the many significant parameters used in mission design studies and data 
reduction analysis is the ray tangent height, which is the minimum altitude above the 
Earth's surface of the ray from the instrument optical center to some point on the 
solar disk, usually the center of the visible disk. The data event time is defined 
as the time interval required for the spacecraft to pass between two specified tan- 
gent heights. For example, on a sunset event, the data collection period may extend 
from an altitude of about 100 km to some small negative altitude. (Refraction being 
neglected, the negative geometric altitude would correspond to a ray just grazing the 
Earth's surface.) 

The duration of the event depends strongly on the orientation of the orbital 
plane of the spacecraft relative to the Earth-Sun line. For example, for a given set 
of tangent point altitudes, the event time will be a minimum when the Sun views the 
orbi t "edge-on • " 

The orbital plane, however, rotates relative to the Sun because of two different 
effects, one dynamic in origin and the other gravitational. If the density of the 
Earth were a function only of the distance from the center of the Earth, the orbital 
plane would remain essentially fixed in inertial space. Thus, as the Earth revolves 
about the Sun in its annual excursion, the orbital plane would appear to rotate rela- 
tive to a solar observer. Because the Earth is not really a perfect sphere gravita- 
tionally, an inertial rotation of the orbit plane results, which is superimposed on 
the annual mean motion. 

One consequence of this motion is that depending on the size of the orbit, the 
inclination of the orbital plane to the Earth's equator, and the declination of the 
Sun, there may come a time when a solar observer would see the orbit just graze the 
edge of the Earth's projected disk. Thus, to an observer on the spacecraft, the Sun 
would neither "rise" nor "set,” but would appear to sink toward the horizon, just 
graze the Earth's surface, and begin to rise again. In fact, depending on the or- 
bital geometry, a condition could arise in which the orbital plane is normal to the 
Earth-Sun line, and hence, for a circular orbit, the Sun would appear stationary to 
an observer on a spacecraft. 

The fundamental parameter on which these conditions are related is the so-called 
"beta angle" (6), which is the angle the Earth-Sun vector makes with its projection 
onto the orbital plane of the spacecraft (fig. 1). 

The following development assumes a general elliptic orbit. However, although 
the principle is exactly the same, much of the algebra becomes complicated very 
rapidly and tends to obscure the basic dependence on 3 * Thus, once the general 
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equations are developed, a restriction to circular orbits will be imposed. The 
simplified equations which result permit a ready physical interpretation which can 
easily be extrapolated to the more complex case of elliptic orbits. 

The basic approach to be taken is to project the radius vector of the spacecraft 
onto a plane normal to the Earth-Sun vector and view the orbit as seen by an observer 
on the Sun. The minimum projected radius of the spacecraft as thus seen, which is 
directly related to B, is then used to determine whether or not event conditions 
exist, and if so, to determine the event times, elevation and azimuth angles, and 
other parameters of use to the mission planner. 

Figure 2 illustrates the geometry. The fundamental plane defined in figure 2 is 
the plane of the spacecraft orbit, with pole w. The vector e Q is a unit vector to 

some part of the solar disk. (The center of the disk is assumed here. ) It is tac- 
itly assumed that the Sun is at an infinite distance. With very little additional 

complexity, the coordinate origin could be transferred to the spacecraft, and a fi- 
nite distance assumed for the Sun. A slight correction to the solar right ascension 
and declination would thereby be introduced to give these quantities as seen from the 
spacecraft coordinate system. These small parallax effects are neglected in the 
present text. If the right ascension and declination of the Sun are given by a 
and 6 g for a specific date and time (see the appendix), then 


cos 6 cos a 
s s 


e s = 


cos 6 sin a 
s s 


sin <5 


( 1 ) 


^ A 

Symbols p and Q represent unit vectors in the orbital plane, with p pointing 
toward periapsis and Q pointing toward a true anomaly of 90°. The symbol w 
represents a unit vector normal to the plane of the orbit along the positive direc- 
tion of the angular momentum vector. In terms of the usual orbital elements (e.g. , 
see Escobal 1965), 


P = 


cos a) cos ft - sin w sin ft cos i 
cos a) sin ft + sin u> cos ft cos i 


sm a) sm i 


( 2 ) 


-sin a) cos ft - cos a) sin ft cos i 


Q = 


-sin u> sin ft + cos u) cos ft cos i 


cos a) sin i 


( 3 ) 
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sin ft sin i 


W = 


-cos ft sin i 


cos 1 


(4) 


where 0 ) is the argument of periapsis, ft is the longitude of the ascending node, 
and i is the inclination of the orbital plane to the Earth equator, 

4 - 

The spacecraft radius vector R is given by (Escobal 1965) 


R = x n P + Y n Q 


(5) 


where 


Xq = a(cos E - e) 

and 


Y 


0 " 



2 . 
e sin 


E 


( 6 ) 


(7) 


In equations (6) and (7), a is the semimajor axis of the orbit, e is the eccen- 
tricity, and E is the eccentric anomaly, related to the true anomaly v and the 
mean anomaly M by the equations 



- e 
+ e 



(8) 


M = n ( t - T) = E - e sin E 


( 9 ) 


The last equation relates E to time t, where T is the last time of periapsis 
passage, and n is the mean angular rate of the spacecraft in its orbit 



In equation (10), y is the gravitational constant of the central planet, which for 
the Earth is 398 600.64 km 3 /sec 2 . 
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Define the orthogonal vector triad 


A A 



( 11 ) 


W X e 

s 

A A 

W X e 

s 


( 12 ) 


(13) 


where e^ is a unit vector to the Sun, and e 2 is a unit vector in the orbital 
plane. 

From the geometry of figure 2, we see that 


sin 3 = W • e g (14) 


cos 3 = W x e s 


(15) 


where 3 is the angle between the Sun vector and the orbital plane. Thus, from 
equation (12), we can write 


A A 

W X e 

A _ s 

e 2 cos 3 

and by use of the vector identity (e.g., Kaplan 1973) 

(A x B) x C = (C • A) B - (C * B ) A 
we can write e^ as 


e sin 3 - W 

A ^ S 

e 3 ~~ cos 3 


(16) 


(17) 
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Let R 2 and R-j denote the projections of R onto e ? and e,, 
respectively. Then 2 3 


R 2 = R * e 2 - 


-r r a. a . 

R . (W x e ) 


cos 6 


and using another vector identity from Kaplan (1973) 


+ . ■> > + 
(A x B) • C = (C x A) • B 


we get 


R 


2 


/V . A 

(R x w) • e 
s 

cos 6 


From figure 2 


A A 


P x w = -Q 


Q x w = P 


and hence, using equation (5), we obtain 


-*0(2 • %) ♦ y 0 (p . g 


R 2 " 


cos 8 


(18) 


Similarly, since R • W = 0 by construction, we get 


R 3 “ [ x 0 ( p * e s ) + Y 0 (2 • e J] tan 6 


(19) 


Define p, the magnitude of the projection of R onto the e 9 -e, plane by 

^ J 


2 2 2 
p = r 2 + r 3 
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which, with the direction cosines of the e s vector 


A A 


p = P • e 


q = Q • e 


w = W • e 


can be written as 


p 2 = R 2 sin 2 B + (X Q q - Y 0 p ) 2 


( 20 ) 


P 2 = R 2 - (X Q p + Y 0 q) 2 


(21 ) 


where 


R 



a ( 1 - e cos E) 


( 22 ) 


and use has been made of the direction cosine property 


2 2 

p + q + 


= 1 


(23) 


Given a particular value of p corresponding to an event altitude (tangent 
point height, for example), equation (20) or (21) can be solved by substitution from 
equations (6) and (7) for E, the eccentric anomaly for which the event occurs* A 
second event condition would give a second value of E, and consequently the time 
interval between the two events can be determined from equation (9) as follows: 


fit 


(E - E ) - e (sin E 2 - sin E ) 
n 


(24) 


Equation (20) or (21) is in general a quartic equation for cos E (or sin E) , 
and hence there are eight roots to be investigated, since the sign of the sine and 
cosine functions only determines one of two quadrants into which E can fall* How- 
ever, only four of the eight values of E will yield the correct value of p, and 
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hence equation (20) must be used to determine which four roots correspond to the 
physically real roots. (See fig. 3 and the numerical example presented later. ) 

Some other parameters of considerable importance in mission design studies are 
the elevation and azimuth angles of the Sun as seen in a local spacecraft coordinate 
frame and the time rates of change of these angles. 

Form two new unit vectors: ^(1) e r = $/|&|, a unit vector along the radius vec- 
tor of the spacecraft, and (2) W x e r (fig. 4). For circular orbits, the velocity 
is parallel to w x e r . Thus, the triad -w, W x e r , and e r form a useful local 
coordinate system, with W x e r pointing "forward," -w pointing out the "right 
wing," and e r pointing "upward" along the positive radius vector, the local 
vertical. (See fig. 5. ) 

In this coordinate system, centered at the spacecraft, any vector, specifically 
e , can be defined by E 0 , the elevation angle above or below the local horizontal 

O A A A Xd 

(the a -W, W x e g plane), and A z , the azimuth angle measured clockwise from, say, 
the W x e r vector in the local horizontal plane. 

From figure 5 and equations (5) and (21), the elevation angle is 



The minus sign is chosen by convention to make the elevation angle negative when it 
is measured below the local horizontal. 

Let a x and a y be the components of e Q in the -w and w x e r directions, 
respectively. Then, 


0 


x “ 


-e c 


A 

• w 


(26a) 


and 


°y = e s • (w X £ r ) = -(e r • e 2 ) cos B (26b) 

From equations (14) and (15), o x and a can be written as 

°x = -sin 6 (27) 

<3*0 - P Y 0 
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Thus the azimuth angle and its quadrant can be determined from 



cr 

y 

cos A„ = — 

/” 2 2 

Fx + °y 

From figure 5, it can be seen that since e^ is a unit vector, 
12 2 

\ o + a = cos E 

V x y Z 

and so equations (29) and (30) appear as 


sin A z 


-sin g 
cos 


cos A z 


q x o - p y o 

R cos E„ 

i 


( 29 ) 


(30) 


(31 ) 


(32) 


(33) 


It is of interest to determine the geographic position of the subtangent point 
at the event time. Let the subscript t refer to either a sunrise or a sunset 
event. The position vector of the spacecraft at the event time is found from equa- 
tions (5) through (7) with the appropriate value of E used in equations (6) and 

(7). Figure (6) shows the tangent point geometry, where p is the tangent point 

^ t 

vector, and is a vector which forms a triangle with R t and p points toward 

the center of the Sun, 

The vector S can be written as 
t a 

S = e R cos 0 
t s t 



1 1 



where 0 is the angle between the radius vector and the Sun vector (see fig* 6) , so 
that 







e 

s 


) 


e 

s 


The right ascension 


and declination 6 ^ 


of p t 


can be found from 


P x = P t cos 6 1 cos a t 


p = p cos 6 sin a 
y t t t 


p = p sin 6 
z t t 




(34) 


(35) 


The geocentric latitude of the subtangent point is identical with the declina- 
tion* To get the longitude, we must compute the sidereal time at the event time. 

The orbital elements used in all the previous calculations are assumed to be 
given at a specified time, t (GMT) and are given with respect to an inertial coor- 
dinate system in which the x-y plane is the Earth equatorial plane, and the x-axis is 
defined by the direction of the vernal equinox* (See Smart 1977, Brooks 1977, or 
Escobal 1965*) The sidereal time, as applied here, locates the Greenwich meridian in 
this coordinate system (see fig* 7) from which we get the relation with reference to 
the sub tan gent point 


a fc = STG + X 


(36) 


where STG is the sidereal time at Greenwich and A t is the local longitude (positive 
east). The Greenwich sidereal time at 0 hours GMT (STGO) on the date is found from 
the equations in the appendix* If the orbital elements are specified at time t 
(GMT) and if At is the time of the event measured from t, then the sidereal time 
of the event is given by 


STG = STGO + 0.25068447 At (37) 


where At is given in minutes from 0 hours GMT, and the constant is the rotational 
rate of the Earth in degrees per mean solar minute. The latitude and longitude of 
the subspacecraft point can be found in a similar manner by using r in place of 
p^ above. t 


Simplification for Circular Orbits 

The imposition of circular orbits permits a great simplification in the evalua- 
tion of the above equations. First, one can set u) = 0, since periapsis is undefined 
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for a circular orbit, and measure E from the ascending node. This permits the 
equations for p and Q to be written as 


P 


cos ft 
sin ft 


Q = 


-sin ft cos i 
cos ft cos i 
sin i 


and equations (6) and (7) reduce to 
Xq = R cos E 


( 38 ) 


(39) 


(40) 


Yq = R sin E 


(41 ) 


where now R is a constant, and E = M. Equations (20) and (21) reduce to 


p 2 = R 2 sin 2 3 + R 2 (q cos E - p sin E) 2 


(42) 


= R 2 [1 - (p cos E + q sin E) 2 ] 


The minimum value of p is found from equation (42) to be 


Pmin = R i sin 6 I 


(43) 


which is the semiminor diameter of the elliptical projection of the circular orbit 
onto the plane. If P m ^ n < R e r the radius of the Earth, there will be a pair 

of distinct sunrise and sunset events. If P m ^ n = R e , the orbit will just graze the 
surface of the Earth as seen by a solar observer, and the sunrise and sunset roots 
coalesce into one value. Finally, for P m ^ n > R e , there will be no sunrise or sunset 
events. 


If we restrict our attention now to the situation which produces a pair of 
events, there will be four physically meaningful roots of equations (42), as stated 
earlier, corresponding to the four points where the spacecraft pierces the cylinder 
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circumscribing the Earth as seen from the Sun and as seen in figure 3. Two of these 
points (roots 1 and 2) correspond to the actual sunrise and sunset conditions, 
respectively, and pierce the "shadow cylinder" of the "shadow side" of the Earth; the 
other two roots (3 and 4) correspond to conditions at the same value of p but 
pierce the shadow cylinder of the "Sun side" of the Earth, These roots can be dis- 
tinguished as follows: R • e g is the component of & lying along the Earth-Sun 

vector. If R • e >0, the spacecraft is on the Sun side of the Earth (i.e., sun- 
ward of the terminator), whereas if £ • e g < 0, the spacecraft is on the shadow side 
(i,e,, behind the terminator). Thus, the real events occur for the two values of E 
for which 


& • e < 0 
s 


(44) 


From equation (5), this expression becomes 
X Q p + Y 0 q < 0 


and for circular orbits 


p cos E + q sin E < 0 


(45) 


The distinction between sunrise and sunset can most readily be made by examining 
the sign of p* From equations (42) 


P = 


-nR 


pq cos 2E + 


2 

q “ 


sin 2Ej 


(46) 


• • 

A positive p value identifies a sunrise event, and a negative p value corresponds 

to a sunset event. 

The values of E which give a specific p can be found by solving equa- 
tions (42) to obtain 


cos E 




t qVi 


2 , 2 .2 

p + q - h 

2 2 
p + q 


(47) 


where, taking into account the condition shown in equation (45), h is the positive 
constant 
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h " + \/’ - (I ) 2 


(48) 


Equation (47) gives two roots for cos E and hence four roots for E. Two of 
the four roots can be eliminated by using equation (42), as they will not produce the 
correct value for p, and the sunrise and sunset roots can be identified from 
equation (46). 

If we write out the expression for sin 3 from equation (14) by using 
equations (1) and (4), 


sin 3 = cos i sin + sin i cos 6 sin (ft - a ) (49) 

s s s 


For a given orbit, 3 is a maximum when the orbit plane is as nearly normal to the 
Sun vector as it can get, that is, when ft - a g = 90° or 27 0°. From equation (49) 
this condition corresponds to 


e 


max 


±(i + « 8 > 


(50) 


For a grazing condition, P m i n ~ R e ' and we have from equation (43) that 


R 

sin e G = ^ 


(51 ) 


where 3q corresponds to the value of 3 for a grazing condition. Thus if 
3 max < the orbit will not experience a grazing situation - there must be events 

for that orbit. Of course, 5 S changes throughout the year from -23.5° to 23.5°, 
and hence there may be grazing possibilities at some times of the year and not at 
others. A grazing condition is most likely to occur near the solstices, and hence 
the maximum inclination an orbit can have in order to insure that it never experi- 
ences a grazing situation can be found from equations (50) and (51) if we set 
6 g = 23.5°. The following table reflects these vaules: 


Orbit altitude, km 

3, deg 

Snax' de 9 

300 

72.76 

49.32 

400 

70.22 

46.98 

500 

68.02 

44.58 

600 

66.07 

42.63 

700 

64.30 

40.86 

800 

62.69 

39.25 

900 

61 .60 

37.76 

1 000 

59.82 

36.38 
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These calculations show that if i < i then no grazing condition can be 

experienced for that orbit, and each orbit will experience a pair of events. How- 
ever, if i > i max * there is no guarantee that a grazing situation will be experi- 
enced, as the grazing condition depends on the angles ft, 6 g , and ot g , all of which 
vary throughout the year. However, since ft and ot g both change more rapid than 
6 , especially near the solstices, the probability is high that a grazing event will 
occur. 


The expressions for the elevation and azimuth angles assume simple forms for 
circular orbits. For these orbits, the unit radius vector becomes 


e = P cos E + Q sin E 
r 


(52) 


and hence equation (25) can be written 


sin E^ = p cos E + q sin E = 



(53) 


Similarly equation (28) reduces to 
0^ = q cos E - p sin E 


Then , using equation (53), 


(54) 



(p cos E + q sin E) 


cos E a 
£ 


(55) 


a result we found earlier from geometric considerations, and the azimuth follows from 
equations (29) and (30) 


sin A z 


-sin 3 
cos E 

£ 


cos A z 


q cos E - p sin E 
cos E 

£ 




(56) 
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The elevation and azimuth angle rates can also be computed from equation (53) 


dE„ 


dt 


n(q cos E - p sin E) 
cos E rt 


= n cos A„ 


(57; 


Comparing equations (55) and (42), we see that 


I 


2 2 p 

a to = — = cos E n 
x y R l 


(58) 


and hence 


sin A„ = 


-R sin 3 


P • 
min 


(59) 


from which 


dA 


-PP . 
mm 


dt 2 

p cos A 


(60) 


Finally, it may be of some interest to compute the rate of change of 3, 
equation (14) 


From 


cos 3 — “ = W • e + W • e 
dt s & 


(61 ) 


If w is the angular rate of the Earth about the Sun, then e can be written 


s s 


-cos e cos 6 sin a - sin e sin 6 
s s s 


cos e cos 6 cos a 
s s 


sin e cos 6 cos 01 
s s 


(62) 


17 



where e is the obliquity of the ecliptic. (See the appendix. ) With equations 
and (4), equation (61) can be put into the form 


cos 6 = -o) [sin i cos e cos 6 cos (ft - a ) 

dt s s s 


+ sin e (sin ft sin i sin 6 - cos i cos 6 cos a )] 

s s s 


+ h sin i cos 6 cos(ft - a ) 
s s 


( 1 ) 


(63) 


The bracketed term in equation (63) results from the annual excursion of the Earth 
around the Sun, and the second term results from the secular perturbation in the line 
of nodes. The angular rate u> g (in degrees per day) can be well approximated by the 
expression 


co s * 1.00288 + 0.03352 cos [0.98562 (d - 2)] (64) 

where d is the day number of the year, and S2 can be found from the perturbation 
equations (Brooks 1977 or Escobal 1965). For circular orbits & can be written as 


h = - 1 J 2 (ir) (cos i)n 

in which J = 1.08228 x 10" 3 is the second-order perturbation constant, and ft has 
the same units as n. 

One final parameter which is of frequent interest might be termed the "skewness" 
of the data profile - this is, the length of the angular arc traced by the subtangent 
point during a sunrise or sunset event. If we let this angle be y, then 


cos y 



(65) 


where the subscripts 1 and 2 refer to the beginning and end of the event, respec- 
tively. Substitute from equation (34), use equation (25), and simplify to obtain 


cos y 


e • e - sin E rt sin E rt 
r 1 r 2 l 2 *1 
cos E cos E 

l 2 £ 1 


( 66 ) 
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Equation (52) permits this to be written as 


cos y 


cos (E^ - E_ 1 ) - sin 


cos 



cos 




( 67 ) 


This arc length, when expressed in units of kilometers on the surface, can range from 

0 km for 6=0° to several hundred kilometers for 8=6 

^max 

Numerical Examples 

Some numerical examples illustrating the computations and some applications 
might clarify the concepts introduced in the text* One complete calculation of a 
pair of events is presented, and then the application of these methods to a typical 
ground- truth experiment is given. 

Calculation of events .- We assume a set of orbital elements presented below: 


a = 6981.2908 km 
e = 0 
i = 57.0° 
a = 266.1083° 
w = 52.5800° 

M = 172.3795° 

Date = November 12, 1985 
T = 0 hours, GMT 


From the equations of the appendix, we find for this date, a g = 227.0949° and 
<$ s = “17.6197°. The sidereal time at Greenwich is 51 .0702°. We assume that sunrise 
or sunset occurs when the tangent ray altitude is -70 km to approximate the effect of 
refraction of the Earth's atmosphere. 

From equations (2)-(4) and (1), we get the unit vectors 


f 0.390311 


P 


-0.63561 


0.66607 
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0.38409 


Q * 


0.76991 
_0. 50962_ 


W 


'-0.83674' 
0.05692 
. 0.54464. 


-0.64885 

-0.69812 

-0,30270 


and the p, q, and w parameters become 


p « -0.01 1 1 4 
q = -0.94097 
w = 0.33832 


Hence, 3 and p . are 
min 

3 = 19.77430° 


^min 


2361.882 km 


Since p min is less than the radius of the Earth (6378 km), there are distinct 
events for this orbit. With p = 6308 km, we get from equations (47), (42), and 
(46), the following table of events for E: 


Root 

E 

p , km 

• 

p, km/sec 

1 

26.407 

6308.0 

-3.002 

2 

333.593 

6372.2 


3 

152.236 

6308.0 

3.002 

4 

207.764 

6240.9 
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and therefore roots 1 and 3 are the physically real roots, with root 1 corresponding 
to sunset (negative p) and root 3 being the sunrise root* The mean angular rate of 
the spacecraft in its orbit is (from eq* (10)) equal to 0.06201 deg/sec. Thus, the 
rise and set times are 


rise 


360 - 172.3795 + 152.2363 
0.06201 


5480.68 sec 


= 1 h 31 m 21 s 


. 360 - 172.3795 + 26.4067 

set ~ 0.06201 


3451 .49 sec 


= 0 h 57 m 31 s 


At sunrise, we find from equations (5)- (7) that 


R = 
t 


'-1 162.106' 
6430.326 
-2457 • 405. 


and hence the tangent point vector is (from eq* (34)) 


p t = 


-3102.970' 
4342.080 
L-3362. 850. 


The right ascension and declination of the subtangent point vector are (from 
eq. (35)) 


a t = 125.551° 


<S t = -32.216° 
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Now, I* 1 31 m 21 s = 91.35 minutes, and thus the Greenwich sidereal time at sunrise is 
(see appendix) 

STG = 51.0702 + (0.25068447)(91.35) = 73.9702° 


and hence the longitude of the subtangent point is (from eq. (36)) 


A t = 125.551 - 73.970 = 51.581° east longitude 


If we repeat this calculation using the sunset root, we get the table below: 


Parameter 

Sunrise 

Sunset 

Time 

Latitude 

Longitude 

1 h 31 m 21 s 
-32.216° 
51 .583° 

0 h 57 m 31 s 
43.910° 
247.614° 


Now, we look at the events at an altitude of 137 km. The range of 137 to -70 km 
(or vice versa) represents a typical altitude range of a Sun-scan-type experiment. 

It is assumed that the Sun's position remains fixed and that the orbital elements 
retain the values they had at time T. By iterating on the results given here, these 
restrictions can obviously be relaxed if higher accuracy is warranted. The calcula- 
tion gives E at sunrise as 156.873° and E at sunset as 21.770°. The sunrise and 
sunset times are I* 1 32 m 35 s and 0* 1 56 m 17 s , respectively. From either pair of times, 
the total event time (from -70 to 137 km, the total time available for measurements) 
is 74 sec. The time spent in Earth shadow is 33 m 49 s . 

To apply these results, the following scheme can be used to determine the range 
of the mission parameters for, say, 1 year. The above calculation is repeated daily 
at 0 hours GMT. The orbital elements are updated each day by using the perturbation 
equations given in Brooks (1977) or Escobal (1965). If this is done for the example 
orbital elements given earlier, then such data parameters as shown in figures 8 
through 11 can be derived. Figures 8, 9, 10, and 11 show as examples of typical out- 
put, the range of subtangent latitude, 8, mission event time, and time in Earth 
shadow, respectively, plotted as functions of days from launch. Figures 10 and 11 
show that mission event time ranges from about 70 to 100 seconds and that time in 
Earth shadow ranges from 29 to 35 minutes except when the orbit approaches a full 
sunlight situation (i.e., near days 60, 130, 240, and 310) and when 6 nears S max 
(66° for this orbit) on days 170 and 350. At these times, the mission event times 
increase dramatically, and the potential for thermal problems in the instrument also 
increases. In addition, the time spent in Earth shadow decreases significantly, and 
there may be a reduction in the instrument recovery time or an adverse effect on the 
data transmission time. 

Calculation of event nearest a specified geographic location. - As a second exam- 
ple of the application of these results to mission planning, suppose that we want to 
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find an event which occurs nearest to some specified geographic location. For exam- 
ple, suppose a balloon launch or ground-based experiment is to be made during the 
sunrise or sunset event which occurs nearest to Laramie, Wyoming, whose coordinates 
are approximately 41° N and 105° W. From figure 8, it is seen that sunrise events 
occur at latitudes of 41° N at about 1 and 21 days after launch, and sunset events 
occur at 30 and 49 days after launch, with corresponding pairs of events occurring 
throughout the year. We want to find the conditions near the Laramie event at day 
21. From a detailed output of the calculations described above, we find the fol- 
lowing table of events for days 21 and 22: 


Parameter 

Day 21 

Day 22 

Date 

Sunrise time 

Subtangent point latitude 
Subtangent point longitude 

November 2, 1 985 
0 h 48 m 06 s 
42.124° N 
87.520° E 

November 3, 1985 
1 h 00 m 27 s 
38.465° N 
83.009° E 


The total time between these events is 24.20583 hours. The orbital period is 
360/0.06201 = 5805 seconds, or 1.6126 hours. Thus there are 24.20583/1.6126, or 15 
events occurring during this interval. Since the longitude change is from east to 
west (opposite to the direction of the Earth's rotation), the total longitude spanned 
at 41° N during this time period is 360 + (87.520 - 83.009), or 364.511°, and thus 
the change in longitude per event is 364.511/15, or 24.300733 degrees per event. The 
change in longitude from 87.520° E to 105° W is about 192°, and thus the event near- 
est to Laramie will be the rounded-off quotient 192/24.300733, or the eighth event on 
November 2. 

During these eight events, the actual longitude change is 194.40586°, so the 
event occurs at 106.886° W, or just west of Laramie. 

The mean time between events is 24.20583/15 = 1.613722 hours, and hence the 
elapsed time for eight events is 12.909776 hours, or 1 2* 1 54 m 35 s . Thus, sunrise 
occurs at 1 2 h 54 m 35 s + 0 11 48 m 06 s , or 1 3 h 42 m 41 s GMT. 

The latitude of the subtangent point is also found by simple interpolation to be 


42.124 


(38.465 - 42.1 24) v 
15 


8 


40.173° N 


A comparison between these interpolated results and a more detailed orbital calcula- 
tion for this specific event gives the following table: 


Parameter 

This case 

Exact 

Time 

Latitude 

Longitude 

1 3 h 42 m 41 s 
40.173° N 
106.886° W 

1 3 h 42 m 40 s 
40.245° N 
106.858° W 
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For most ground-truth experiments tied to satellite rise or set events, the interpo- 
lated results appear to yield results of acceptable accuracy. 


CONCLUDING REMARKS 

An analytical method is developed for determining the geometrical parameters 
which are needed to describe the viewing angles of the Sun relative to an orbiting 
spacecraft when the Sun rises or sets with respect to the spacecraft. These equa- 
tions are rigorous and are frequently used for parametric studies relative to mission 
planning and for determining instrument parameters. 

The text is wholly self-contained in that no external reference to ephemerides 
or other astronomical tables is needed. Equations are presented which allow the com- 
putation of Greenwich sidereal time and right ascension and declination of the Sun 
generally to within a few seconds of arc, or a few tenths of a second in time. 


NASA Langley Research Center 
Hampton VA 23665-5225 
June 11, 1986 
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APPENDIX 


COMPUTATION OF SOLAR POSITION 

The following algorithm will produce the right ascension and declination of the 
Sun to within a few seconds of arc, as evidenced from a considerable amount of com- 
parison with solar position data published in various annual almanacs. This appendix 
is completely self-contained. No other tables or reference materials such as alma- 
nacs or ephemerides are needed to compute the parameters discussed in the text; 
hence, these relations are quite suitable for computer use. Angles are computed mod 
(360). The steps and equations of the algorithm are 

1. Compute the Julian date (Almanac for Computers 1 980) 


JD = 367K - + <( - M - 4 + 9 >/ 12> ] ^ + + i + 1721013.5 + 

where 

JD Julian date 

K year (e.g., 1984) 

M month (1 < M < 12) 

I day of month (1 < I < 31 ) 

GMT Greenwich mean time, hours (0 < GMT < 23.99) 

The symbol < > denotes the maximum integer value. 

2. Compute the sidereal time at Greenwich at 0 hours GMT of a given date as follows 
(Escobal 1965): 

STGO = 99.6909833° + 36000. 7689°T + 0.00038708°T 2 

u u 


where T u is defined by 


JD - 2 415 020 
u 36 525 


The equation for STGO includes the precession of the equinoxes but neglects the nuta- 
tion terms, that is, the equation of the equinoxes. Hence, this is the mean sidereal 
time at 0 hours GMT. For any GMT, the sidereal time is then found from 


STG = STGO + 0. 25068447GMT 
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where GMT is in minutes 


4. 


5. 


6 . 


Compute the right ascension of the mean Sun (RAMS) (Escobal 1965) 

RAMS = 279.6966778° + 36000. 76892°T + 0.0003025°T 2 

u u 

Compute the obliquity of the ecliptic e (Escobal 1968) 

e = 23.45229444° - 0.0130125°T - 0. 000001 6389°T 2 

u u 

Compute the Sun’s apparent mean anomaly M g (Escobal 1968) 

M = 358.475844° + 35999. 04975°T - 0.00015°T 2 - 0.0000033333°T 3 

S U u U 

Compute the eccentricity of the Earth's orbit e 0 (Escobal 1968) 

e = 0.01675104 - 0.00004180T - 0.0000001 26T 2 

e u u 

Compute the equation of time (Smart 1977) 

E m = y sin 21 - 2e_ sin M + 4e y sin cos 2 £ 

T 2 e s e J s 

1 2 5 2 f 5 2 

- — y sin 4£ — 7 e sin 2M + — ye sin 2M cos 21 

2 2 4e s L 2 ^ e s 

- 4y 2 e^ sin cos 4£ + 77 y 3 sin 6J l + 4“ e 3 sin M 

J e s 3 4e s 


1 3 
12 


sin 3M 


where 


i = RAMS 
y - ta„2 | 


The term in brackets is of third order in the Earth eccentricity and may be omitted 
unless extreme accuracy is required. 
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If one substitutes the values of y and e 0 for mid- 1985, the following equa- 
tion (Smart 1977) will give reasonably accurate numerical values for the decade of 
the 1980's: 


E t = -0.4329 sin £ - 1.7900 cos £ + 2.4848 sin 2£ 

- 0.0083 cos 2£ + 0.0179 sin 3£ + 0.0804 cos 3£ 

- 0.0529 sin 4£ + . . . 

The units are degrees. 

8. Compute the right ascension of the true Sun a g (Smart 1977) 
CX S - RAMS - E t 


9. Compute the declination of the Sun 6 g 


tan 6 S = sin a g tan e 


10. If needed, compute the hour angle of the mean Sun (HAMS) and the hour angle of 
the true Sun (HAS) (Smart 1977) 


HAMS = GMT - 1 2 h - Longitude 
HAS = HAMS + E t 


Longitude is defined as positive east of Greenwich. The west longitude of the 
true Sun is just the hour angle of the true Sun measured at Greenwich. The latitude 
of the true Sun is, of course, identical with the declination. 
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Figure 1 Definition of the angle 0. W is a unit 
vector normal to the orbital plane. 









Figure 7.- Relationship between Greenwich meridian, longitude, 
and sidereal time at Greenwich, STP is the subtangent 
point. 
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